Coulomb blockade double-dot Aharonov-Bohm interferometer: harmonic 
decomposition of the interference pattern 
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For the solid state double-dot interferometer, the phase shifted interference pattern induced by 
the interplay of inter-dot Coulomb correlation and multiple reflections is analyzed by harmonic 
decomposition. Unexpected result is uncovered, and is discussed in connection with the which-path 
detection and electron loss. 
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Introduction. — As an analogue of Young's double-slit 
interference [l|, a ring-like Aharonov-Bohm (AB) inter- 
ferometer with a quantum dot in one of the interfering 
paths is of interest for many fundamental reasons and 
receives extensive studies 0, II H- Ver y recently, an 
elegant study was further carried out for the closed-loop 
setup, with particular focus on the multiple-reflection 
induced inefficient which-path information by a nearby 
charge detector @. 

In this report we consider an alternative solid-state 
AB interferometer, say, electron transport through par- 
allel double dots (DD) in Coulomb blockade regime, 
as schematically shown in Fig. 1. Existing studies 
on this DD setup include the cotunneling interference 
@, H, H, H3] j and two-loops (two fluxes) interference with 
the two dots as an artificial molecule [TJ ^ n this 

work we restrict to the DD without tunnel-coupling be- 
tween the dots, which enables us to focus on the interplay 
of Coulomb correlation and quantum interference. 

At first sight, this electronic setup looks similar to the 
double-slit interferometer [l|. Nevertheless, the underly- 
ing solid-state specifics will result in such novel behav- 
ior as non-analytic current switching tuned by very weak 
magnetic field [131 ] . In present work, we extend our previ- 
ous study 13] to more general situations for off-resonant 
DD levels, in the presence of which-path detection and 
lossy channels. We find the interference pattern to be 
asymmetric with respect to the magnetic field, and at- 
tribute it to the underlying Coulomb correlations. More- 
over, similar to Ref.0, we present an analysis of harmonic 
decomposition for the interference pattern. We find that 
the phase shifts of the harmonic components of current 
are linearly scaled with the order of harmonics, and the 
high-order harmonics, in contrast to Ref. |f| will be de- 
phased by the which-path detection. 

Model and Method. — Consider the double dots con- 
nected in parallel to two leads. For simplicity we assume 
that in each dot there is only one level, £1(2) , involved 
in the transport. Also, we neglect the spin degrees of 
freedom. In the case of strong Coulomb blockade, the 
effect of spin can be easily restored by doubling the tun- 



FIG. 1: Schematic setup of a double-dot Aharonov-Bohm in- 
terferometer. The solid-line trajectory is the usual first-order 
harmonic interference, while the dashed-line trajectory repre- 
sents contribution to higher(second)-order harmonic interfer- 
ence. To address dephasing and electron-loss effects, a nearby 
quantum-point-contact (QPC) detector and lossy channels to 
side-reservoirs (with rates 71(2)) are introduced. 



neling rates of each QD with the left lead. The system is 
described by the following Hamiltonian 

H = H + H T + E » d l d >i + Ud\d x d\d 2 . (1) 

Here the first term, H = Y,k[ E ^La\ L a kL + E kR a\ R a kR ], 
describes the leads and Ht describes their coupling to the 
dots, 



(2) 



fi .k 



where fi = 1 , 2 and a kL and a kR are the creation opera- 
tors for the electrons in the leads while d\ 2 are the cre- 
ation operators for the DD. The last term in Eq. (Q]) de- 
scribes the interdot repulsion. We assume that there is no 
tunnel coupling between the dots and that the couplings 

are independent of energy. 



of the dots to the leads, t 
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In the absence of a magnetic field one can always choose 
the gauge in such a way that all couplings are real. In 
the presence of a magnetic flux <I>, however, the tunneling 
amplitudes between the dots and the leads are in general 
complex. We write = ^(ffje^W , where 

is the coupling without the magnetic field. The phases 
are constrained to satisfy 4>\l + — 4>2L — 4>2R = 4>, 
where </> = 27r$/$ - 

To account for dephasing effect, we introduce a which- 
path detection by a nearby quantum point contact 
(QPC)[1], with a model description as in Ref. 0. To 
make contact with conventional double-slit interferome- 
ter, we also introduce electron lossy channels. Slightly 
differing from Ref. [1, instead of the semi-infinite tight 
binding chain introduced there, we model the lossy 
channels by attaching each dot with an electronic side- 
reservoir, which is particularly suited in the master equa- 
tion approach. The side-reservoir model was originally 
proposed by Biittiker in dealing with phase-breaking ef- 
fect [IH, i.e., electron would lose phase information after 
entering the reservoir first, then returning back from it. 
But here, we assume that the reservoir's Fermi level is 
much lower than the dot energy. As a result, electron 
only enters the reservoir unidirectionally, never coming 
back. 



Here p^ denotes the time derivative of the "n"- 
conditional DD state. Owing to the neglected spin de- 
grees of freedom in constructing the Hilbcrt space, as a 
compensation, we have replaced with 2r^ to equiv- 
alently restore the spin effect. In the above equations, 
A = Ei — E 2 is the disalignment of the DD levels. 
t l(r) = 2nD L(R) \t L(R) \ 2 , and j 1{2) = 2ttD 1(2) \t 1{2) | 2 , 
are the respective coupling rates of the DD to the left 
and right leads, as well as to the side reservoirs. Dl(r) 
and -Di(2) are the density of states of the leads and reser- 
voirs, while i_L(_R) and ti( 2 ) are the respective tunneling 



The transport properties of the above described sys- 
tem can be conveniently studied by the number-resolved 
master equation [H, El El- The central quantity of 
this approach is the number-conditioned reduced state, 
p^ (t) of the double dots, where n is the electron num- 
ber passed through the junction between the DD and an 
assigned lead where number counting is performed. Very 
usefully, p( n '(t) is related to the electron-number distri- 
bution function, in terms of P(n,t) = Tr[p( n \t)], where 
the trace is over the DD states. From P(n, t) the cur- 
rent and its fluctuations can be readily obtained. For 
current, for instance, it simply reads I(t) = ed(n(t))/dt, 
where (n(t)) = ~^2 n nP{n,t). In practice, instead of di- 
rectly solving P(n,t), much simpler equation-of-motion 
technique is available for the calculation of current and 
current fluctuations [Tt], Gil • 

In large bias sequential tunnelling regime and under 
inter-dot Coulomb blockade (i.e. the DD can be occupied 
at most by one electron) , the Hilbert space of the DD is 
reduced to |0) = |00), |1) = |10), and |2) = |01), where 
1 10) means the upper dot occupied and the lower dot 
unoccupied, and other states have similar interpretations. 
Following Ref. [H, the "n" -resolved master equation in 
this basis can be straightforwardly carried out as 



(3a) 



(3b) 



(3c) 



(3d) 



(3e) 

I 

amplitudes. For simplicity, in this work we assume that 
7i =72 =7, and Tl = T R = T in the following nu- 
merical results. Finally, 7^ is the dephasing rate between 
the two dots, caused by the which-path measurement of 
QPC [3. 

Phase-Locking Breaking. — In the absence of electron 
loss, i.e., 7 = 0, simple expression for the steady-state 
current is extractable: 

r 2( 7d + 2r R )(l-cos0)-4Asin0 1 } ~* 

1 7 d(7d + 2r R ) + 4A2 + / J ' U 



p® = -2T ip W + ( 7 + T R ) P [r 1] + (7 + T«) pf 1} + e^-^TnpiT^ + e^-^vj^ 
= TlP^ - (7 + Tr) p$> - \e^-^T R p^ - l -e^-MY R p^ 

P2H = ^ - (7 + T R ) - ie«*~-*«>r B p<2> - l -e^-^T R p%> 
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FIG. 2: (A) Phase shifted Aharonov-Bohm interference pat- 
tern, as a consequence of the interplay of Coulomb correlation 
and quantum coherence. Inset: shift of current peak versus 
the level detuning A. Effect of dephasing between the two 
dots is demonstrated by varying 7^. (B) Effect of electron loss. 
With increasing the lossy strength 7, conventional double-slit 
interference pattern is eventually restored. 




FIG. 3: Harmonic decomposition of the interference current 
I(4>)- (A) and (B): Amplitude T n and phase shift j3 n of the 
n-th harmonic component, for different energy detuning A. 
In (B) an excellent linear-scaling relation for the phase shift, 
say, /3„ = n(3i, is observed. (C) and (D): Effect of which-path 
detection on individual harmonic components. In contrast to 
Ref. [H, here the higher order harmonic components are also 
dephased (even more seriously). Interestingly, from (D) we 
find that the scaling relation /3 n — n/3i still holds even in the 
presence of depahsing. 



where Jo = 4rLrR./(4rL + Tr), is the current in the ab- 
sence of magnetic flux. However, in the following we will 
use the current of transport through a Coulomb-blockade 
single dot, h = 2r L r R /(2r L + T R ), to scale the double- 
dot current, in order to highlight the interference fea- 
tures. 

In the absence of dephasing, 7^=0, from Eq. (|4|) we 
have / = / o A 2 /{A 2 + /o[rfl(l-cos0)-Asm(/)]}. Then a 
remarkable switching effect follows this result: as A — * 0, 
I = Iq for <f> = 2wn, while 1 = for any deviation of <f) 
from these values. Detailed interpretation for this novel 
behavior using a SU (2) transformation is referred to Ref. 

m 

If A 0, the Coulomb correlation is manifested in an 
alternative way, see Fig. 2(A), where the current peaks 
deviate from <j> = 27m. For small A, we find that the cur- 
rent is peaking at cf> » A/Tr + 27m, and with a magni- 
tude J max s» (I /la — 1/2Tr)~ 1 . Increasing A, the current 
peak moves rightward. This kind of peak shift implies 
nothing but the symmetry breaking with respect to the 
magnetic-field inversion. This behavior contrasts with 
the following statement (l9j : based on current conserva- 
tion and time-reversal invariance, the Onsager relation, 
say, the symmetry relation of transport coefficients under 
inversion of magnetic field, will lock the current peaks at 
<j> = 27m, for any two-terminal linear transport. This is 
usually referred to as phase locking. Beyond linear re- 
sponse, however, phase locking does not necessarily hold 
in general. 



In Ref. [t| for instance, for the AB setup with a single 
dot embedded in one of the arms, it was found that the 
phase locking is indeed broken under finite bias voltage 
and only in the presence of electron-electron interaction. 
However, for the similar interacting DD system as con- 
sidered here, breaking of phase locking was not found 
under finite bias voltage, in the case of A = and within 
the cotunneling transport mechanism Q. Here, we find 
that the current is asymmetric only for A^O, while it is 
still symmetric for A = 0, being somewhat in agreement 
with Ref. [9|, despite the different transport mechanism 
and dephasing involved here. 

Dephasing and Lossy Effects. — In Fig. 2(A), we also 
observe that the current at 4> = 2n7r is not affected by 
A and 7^. This feature implies that the current in case 
of complete dephasing is the same as the one from con- 
structive interference (i.e. with <fi = 27m), being different 
from the conventional double-slit interference [1]. Wc 
attribute this feature to the closed boundary condition, 
under which the multiple reflection plays essential role. 
By introducing lossy channels, i.e., allowing electron loss 
from the DD to surrounding environment, in Fig. 2(B) we 
see that all these features disappear and the conventional 
double-slit interference pattern is restored by increasing 
the lossy rate 7. The reason is that, as the dots become 
more and more open, the side reservoirs would reduce 
the occupation probability on the dots, thus make the 
Coulomb correlation and back- reflection less important. 

Harmonic Decomposition. — To highlight the effect 
of the closed nature and Coulomb correlation, we fur- 
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ther expand Eq. (4) into Fourier series, I(<f>) = Zq + 
S^Li -^n cos(n</> + /3 n ). Physically, the n-th order har- 
monic stems from trajectories which have n-turns dif- 
ference surrounding the magnetic flux between the two 
interfering partial waves (see Fig. 1). In Fig. 3 (A) and 
(B), we plot the amplitude X n and phase shift n for co- 
herent case. The higher order harmonics are caused by 
the multiple reflections under closed boundary condition, 
while the phase shift is related to Coulomb correlation be- 
tween the two dots (note that there is no phase shift for 
nonintcracting DD). Quite surprisingly, in Fig. 3(B) we 
observe that (3 n = n(3\. At this stage, unfortunately, we 
failed to find a satisfied interpretation for this interesting 
result, so would remain it for future study. 

In Fig. 3 (C) and (D), dephasing effect is shown. Note 
that, as explained in the model description, in this work 
dephasing is modelled by a which-path detection. Inter- 
estingly, we observe here that the which-path detection 
also dephases the higher order harmonic components of 
current, unlike in Ref. [f| where only the first-harmonic 
amplitude is strongly reduced by the detection, while the 
second one is almost insensitive to such detection. We 
understand this discrepancy as follows. Taking the sec- 
ond harmonic as an example, the first partial wave has 
an amplitude oc t\e %xi , while the second partial wave 
oc t2e lX2 ■ t\e tXl ■ t2e %X2 . Here, we formally denote the 
transmission amplitude through dot-l(2) by ti(2)> which 
is in general complex, e.g., containing the Aharanov- 
Bohm phase; Xi(2) an d Xi(2) are the random phases 
caused by the charge detection when the electron passes 
through dot-l(2) in the first- and second-order trajecto- 
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